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Abstract. We present a comprehensive study of the observational dependence of the mass-loss rate in stationary stellar winds 
of hot massive stars on the metal content of their atmospheres. The metal content of stars in the Magellanic Clouds is discussed, 
and a critical assessment is given of state-of-the-art mass-loss determinations of OB stars in these two satellite systems and the 
Milky-Way. Assuming a power-law dependence of mass loss on metal content, M oc Z'", and adopting a theoretical relation 
between the terminal flow velocity and metal content, Voo Z" '^ (Leitherer et al. 1992), we find m = 0.83 ± 0.16 for non- 
clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than lO'^Lg. 
Within the errors, this result is in agreement with the prediction m = 0.69 ± 0.10 by Vink et al. (2001). Absolute empirical 
values for the mass loss, based on Ho- and ultraviolet (UV) wind lines, are found to be a factor of two higher than predictions in 
this high luminosity regime. If this difference is attributed to inhomogeneities in the wind, and this clumping does not impact 
the predictions, this would imply that luminous O and early-B stars have clumping factors in their Ha and UV line forming 
regions of about a factor of four. For lower luminosity stars, the winds are so weak that their strengths can generally no longer 
be derived from optical spectral lines (essentially Ho-) and one must currently rely on the analysis of UV lines. We confirm that 
in this low-luminosity domain the observed Galactic WLR is found to be much steeper than expected from theory (although 
the specific sample is rather small), leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 
at luminosities of L ~ 10'* '^Lq, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot 
luminous stars would turn out to be overestimated as a result of wind clumping, but the degree of clumping would be rather 
independent of metallicity, the scalings derived in this study are expected to remain correct. 
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1. Introduction 

Understanding the properties and evolution of massive stars 
in low metallicity environments is of fundamental importance 
in astrophysics. This is amply illustrated by referring to the 
anticipated role of massive stars in the early universe, where 
they are thought to be responsible for its re-ionization (e.g. 



[Haehnelt et ar][200T] [Wyithe & L oeb"2003'), galaxy forma- 
tion, and chemical evolution of young galaxies and of the intra- 
galaxy medium. Strong indications of gamma-ray bursts occur- 
ring primarily at low metallicity (e.g. [Gorosabel et al.||2005[ l 
and massive stars being progenitors (e.g. Hjorth et al. |2003 1, 
adds to the point. To a large extent, the evolution of massive 
stars is controlled by the effects of mass loss, therefore mass 
loss has a direct and/or indirect impact on all the phenomena 
mentioned above. This makes a quantitative handle on the de- 
pendence of mass loss on the metal content of the environment 
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out of which these early stars form so vital. Obviously this ar- 
gument also applies to massive stars in the local universe. 

In the last decade the development of realistic stellar at- 
mosphere models has allowed accurate determination of the 
wind parameters of a sizeable sample of early-type Galactic 
stars ( e.g. [Puis et al.|1996[|Herrero et al.|2002[[Repolust et al. 



Table 1. Overview of the ions that dominate the line driving in 
early -type stars - for solar metallicity with - 10^ Lq. See 
Sect.|2]for a discussion. 



|2004[ ICrowther et al.||2006| l. This, in turn, has allowed us to 
test our understanding of aspects of stellar wind hydrodynam- 
ics, including its driving mechanism, to an unprecedented level. 
During this period new instrumentation developements on large 
telescopes have resulted in significant numbers of extragalac- 
tic stars being spectroscopically studied. Examples of extensive 
surveys are those of Massey et al. (2004, 2005) and of Evans 
et al. (2005), although there are many others and the relevant 
ones will be discussed later. The latter program is the VLT- 
FLAMES Survey of Massive Stars, in which close to 100 O- and 
early B-type stars were observed in the Magellanic Clouds. The 
inclusion of line blocking/blanketing in the modeling of hot 
star atmospheres (e.g. |Hubeny & Lanz|1995t|Himer & Miller| 
1998| l allows for a consistent description of the effects from 



vast numbers of overlapping metal lines, particularly regard- 
ing a detailed treatment of the effects of chemical composition, 
and has paved the way for a robust quantitative comparison of 
the wind strengths of hot massive stars in different metallic- 
ity environments. The most recent step forward in the model- 
ing technique has been the development of an automated fitting 
method, opening up a means to analyse large samples in a ho- 
mogeneous way ( (Mokiem et al.|2005 1 



This paper provides a critical assessment of the current 
standing of mass-loss determinations in the Galaxy and the 
Large and Small Magellanic Clouds, and is specifically aimed 
at establishing the dependence of mass loss on the chemi- 
cal abundance pattern, notably metal content Z. The Small 
Magellanic Cloud (SMC) with Z at about l/5th of the Galactic 
value is the lowest metallicity we can so far probe in reasonable 
detail. The Large Magellanic Cloud (LMC) with approximately 
half the metal content of the Galaxy provides an intermediate 
environment. 

The only possibility, in the forseeable future, of studying 
stars in environments with significantly lower metallicity than 
the SMC is the massive stellar populations of the Local Group 
dwarfs GR8, Sextans A and Leo A. These three galaxies have 
metallicities betweem 0.1 and 0.05 Z©, but their low starfor- 
mation rates and distances of around 1-2 Mpc makes studies of 
large numbers of young massive stars extremely difficult. Given 
this difficulty of establishing the dependence of mass-loss rates 
upon metallicity below 1/5 Zg, we must rely on predictions to 
access this interesting part of parameter space. Although we 
do not deal with the theory of wind driving mechanism, we 
will compare predicted and observed wind strengths to estab- 
lish successes and failures of the theory such that we may iden- 
tify aspects of the M(Z) problem that require further study, and 
regimes in parameter space that seem sufficiently under control 
to allow us to venture out to lower metallicity regimes and the 
early universe (see also |Vink & de Koter^2005[ ). 

In Sect. |2] we briefly review the mass-loss mechanism 
of early-type stars. The present day chemical composition of 
the Galaxy and Magellanic Clouds is discussed in Sect. |3] 
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Sections]?] and |5] describe the methodology used to determine 
mass-loss rates and how to compare these for different chem- 
ical environments. The observed mass-loss relations are pre- 
sented in Sect. ]6] On these relations the global metallicity de- 
pendence of M is determined in Sect. ]7] Finally, in Sect. Js] we 
discuss the observed mass-loss relations in terms of the so- 
called weak wind problem and wind clumping, and end with 
concluding remarks. 



2. The mass-loss mechanism of early-type stars 

We briefly review those aspects of the wind driving mechanism 
of massive early-type stars that are relevant in the context of es- 
tablishing an empirical relation between mass loss and chemi- 
cal composition. For a more in-depth treatment of the physics 
of mass loss, see e.g. Kudritzki & Puis ( 2000| l and Vink et al. 
(2001). The basic mechanism driving the winds of hot mas- 



sive stars is the transfer of momentum from photons to the 
atmospheric gas by line interactions. The driving mechanism 
implies that the properties of the stellar wind will depend on 
the number of photons per second streaming through the pho- 
tospheric layers (reflecting the stellar luminosity), and on the 
number and ability of lines being available - in particular at 
wavelengths around the photospheric flux maximum - to ab- 
sorb or scatter these photons. 

The dependence of the wind driving on the number of lines 
present suggests that mass loss is a function of elemental abun- 
dance. Whether this is indeed the case formally depends on the 
nature of the driving lines. In the hypothetical case that one 
would be in a regime of abundances for which all lines effec- 
tively contributing to the line force are optically thick, mass 
loss would not be a function of elemental abundance. This 
regime is not encountered in even the most metal rich environ- 
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ments known ( Vink et al.|2001 1. In reality, it is found that the 
lines driving the wind are a mixture of optically thin and opti- 
cally thick lines. Representing the distribution of line strengths 
by a power law, one predicts, for Galactic O stars, a ratio of 
the line acceleration from optically thick lines to the total line 
acceleration of a ~ 2/3 (see [Puis et al.|[2000) l, where an even 
larger contribution by optically thick lines would increase this 
value (and vice versa). This ensures a dependence of mass loss 
on elemental abundance. 

Which elements dominate the line force? The answer de- 
pends on the effective temperature of the star (More precisely: 
on the radiation and electron temperature in the wind accel- 



eration regime; see also Vink et al. 1999 Puis et al. 2000| 



Although hydrogen and helium are by far the most abundant 
elements, their impact on the wind driving is modest. Decisive 
for whether or not a species is a significant contributor to the 
line driving is the product: abundance x ionisation fraction X 
number of effective lines. Very roughly, the elemental abun- 
dance times the ionisation fraction of hydrogen and helium are 
similar to that of metals that are in the dominant stage of ionisa- 
tion. It is therefore the number of driving lines that is decisive. 
As H and He - due to their simple atomic structure - have only 
few lines that can effectively contribute to the line force, it is 
relatively abundant complex atoms that are the main contribu- 
tors. 

We estimate the relative contributions of the different ele- 



ments by the use of a Monte Carlo method (Abbott & Lucy 



[1985) |de Koter et aT|[T997) |Vink et al.|[T999; that calculates 
the total momentum transfer from the radiation field to the out- 
flowing gas particles. Although the relative contribution to the 
line force is depth-dependent (see below), here we simply reg- 
ister the atomic number of the elements with which the photons 
interact somewhere in the wind, and we present the results in 
Tab.[T] For a late-O dwarf of solar composition, CNO accounts 
for some 15 percent; iron contributes some 25 percent. Other 
iron-group elements (for instance Cr, Mn, Co, Ni) add a few 
percent, a-elements, such as Ne, Mg, Si, S, Ar, and Ca account 
for about 30 %, with Si being the main contributor with ~ 20 %. 
We mention a decomposition of the line force in terms of iron 
group and a elements as their nucleosynthetic origin is differ- 
ent. The former are mostly produced in thermonuclear Type 
la supernovae, the latter predominantly in core-collapse super- 
novae of types II and Ib/c. Beware that due to their different 
electronic structure specific groups of elements have a differ- 
ent line-strength statistics, and dominate the line acceleration 
at different depths: iron group elements have a somewhat larger 
influence on the mass-loss rate, whereas lighter ions dominate 
the acceleration in the outer wind, thus controlling the terminal 



velocity (cf. Vink et al. 1999^ Puis et al. 2000|. These more sub- 
tle effects are not reflected in the statistics provided in Table [T| 
Mixing of CNO-cycled material to the surface during the 
supergiant phase may affect the relative abundances of these 
three elements. However, in terms of their contribution to the 
line force not much will change when this happens, as these 
three elements have more or less equal numbers of effec- 
tive driving lines near the photospheric flux maximum and 
the C-I-Nh-O abundance remains unaffected by the CNO-cycle 
( IVink & de Koter,2002| ). 



Accounting for the fact that we are interested in the gross 
dependence of wind parameters on metallicity, and particularly 
in the product of mass-loss rate and terminal velocity (see be- 
low), we conclude that a straight mean ofC+N+0, a-elements, 
and iron-group elements is an appropriate abundance value to 
use as aZ measure for O stars. For stars of spectral type B and 
A the contribution of iron increases to up to 50 percent. 

3. Metallicity determinations of early-type stars 

Having identified the elements that dominate the line force, we 
now address the present-day chemical composition of hot mas- 
sive stars in the Galaxy and the Magellanic Clouds. 

Spectroscopy of H ii regions and supernova remnants have 
been used to probe the gas phase composition of the Magellanic 
Clouds (for reviews see for instance Pagel et aL[l978^ |Dufour| 
19901 [Russell & D6pita|[T990l [Garnett[|1999| l. The results de- 
rived from these studies do not necessarily represent the ini- 
tial chemical composition of the stars in these environments 
because of fractionation of some elements from the gas phase 
onto solid state particles (or "dust"). Also, emission line stud- 
ies may be affected by spatial inhomogeneities in the nebula. 
Therefore, to obtain the appropriate chemical composition of 
the winds we are studying it would be pertinent to determine 
the photospheric abundances of the massive stellar populations 
in the Clouds. As the comparison between theoretical and ob- 
served wind strength is usually not done on an individual ba- 



sis (but see |de Koter et al.[[T997 1 it is relevant to specify how 
the observed metal abundance Z - so far the abundance input 
parameter for mass-loss predictions - is derived. If Z is de- 
rived from CNO abundances in either evolved or rapidly ro- 
tating objects one should be careful. In supergiants and even 
giants ( Kom et al.|2000 Lennon et al.|2003" i rotational mixing 
and ejection of the outermost envelope by mass loss may bring 
CNO equilibrium material to the surface. In stars rotating at 
about half of break-up at the surface or more, abundance alter- 
ations due to rotational mixing may already occur very early 
on in the stars evolution ( Yoon & Langer|2005 1. 

The optical spectra of O stars show few spectral lines of 
heavy elements. B stars are relatively rich in absorption fea- 
tures due to C, N, O, Mg, Al, Si, S, and Fe providing a sound 
basis for establishing the metallicity, although we note that the 
Fe lines are intrinsically quite weak (especially so in the SMC; 
Rolleston et al.|2003) . The advantage of B and later- type dwarf 
and giant stars is also that they may be analysed accurately us- 
ing non-LTE line-blanketed hydrostatic atmospheres, without 
winds. 

A major goal of the VLT-FLAMES survey of massive stars 
(Evans et al. 2005) was to observe a large number of such B- 
type stars in the same vicinity as the O-stars to determine their 



photospheric abundances. Hunter et al. (20071 and Trundle 



[et al. ( 2007 1 have presented the chemical compositions of 
34 and 53 B-type stars in the SMC and LMC respectively. 
These are the narrow lined stars from the FLAMES survey, 
with the highest signal-to-noise spectra, analysed with the line- 
blanketed non-LTE code tlusty. These authors determined the 
best estimates for the chemical composition of the Clouds from 
a comparison of results from stellar and nebular work, and we 
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Table 2. Present-day chemical composition of the LMC and SMC from B-type stars. The values are taken from Hunter et al. 
( 2007| l and Trundle et al. (2007) apart from Al and S which are LTE results from Rolleston et al. (2002 2003) 1. For comparison 
the solar abundances of |Asplund et al.| ( [2005| l also given. 
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summarise and extend their results in Table|2] The non-LTE 
absolute abundances for C, N, O, Mg, and Si are quoted, along 
with the simple difference between these absolute values and 
the current standard solar abundances of Asplund et al. ( 2005[ l. 
However detailed non-LTE model atoms for use in calculating 
line profiles for Al, S and Fe in tlusty are required to allow 
reliable line abundances to be calculated consistently in non- 
LTE. The results for Fe are taken from [Trundle et al!] ( |2007| l in 
which a non-LTE atmosphere and LTE line formation calcula- 
tions were used. The results for Al and S are taken from the 
LTE B-star analysis of Rolleston et al.| ([2003 ) and Rolleston 
et al. (20021. Because we are less certain of the validity of 
the absolute values of the LTE results, we quote the differen- 
tial abundances from these papers which were determined from 
analysing Galactic B-stars with identical atmospheric parame- 
ters. We should also note that the [Rolleston et al!] ( |2003] l abun- 
dances come from the analysis of a single star AV304, and the 
LMC results are from 5 stars with spectra of modest signal-to- 
noise. 

The SMC abundance of the a elements (O, Mg, Si) has a 
mean relative to solar of A[q'/H] = -0.7 + 0.1 dex. The dif- 
ferential LTE result of Al is in good agreement with this value, 
although the one S abundance from the star AV304 is some- 
what higher. The Fe abundance is A[Fe/H] - -0.57 ±0.16 dex, 
which, within the uncertainties, is in agreement with the a 
elements. The SMC Fe abundance has also been determined 
through studies of AFGK supergiants, which have more and 
stronger metal lines in their spectra and for which non-LTE 
effects are less important. They tend to show very good agree- 
ment at A[Fe/H] = -0.60+0.1 dex ( |Venn|1999[ ). Differences in 
modelling techniques/codes, diagnostic lines, and atomic data 
will clearly lead to systematic differences between the sets of 
results. The discrepancies between different analyses are be- 
coming smaller, but have not completely been resolved and 
progress is required in the model atmosphere and line forma- 
tion codes, particularly with regard to atomic data for S, Al 

and Fe for B-star analysis. We adopt A[Z/H] 0.7 ± 0.1 for 

the SMC, but will discuss the effects of a potentially different 
metallicity later (see Sect.[7.2|l. 



A[a/H] - -0.36 ±0.1 dex. An extensive comparison of abun 
dance determinations in the LMC was compiled by [Rolleston 



et al. (2002). Reasonable agreement is found between six in- 
dependent studies of the overall metallicity of the LMC - irre- 
spective of the class of object used to trace the chemical com- 
position or the spatial location of the investigated stars. The 
mean metalUcity, based on O and Si, is (A[0/H] + A[Si/H)/2 = 
-0.30 ± 0.08 dex. Some previous analyses using OB stars 
seem to indicate that iron is less deficient than the a-elements. 
Rolleston et al. ( 2002[ l report an a to iron ratio A[a/Fe] ~ 
-0.14 dex for main sequence OB stars, consistent with that ob- 
tained for evolved B-stars ( |Korn et al.|2000| ). However |Trundle[ 
et al.[ ( [2007| l determine a mean differential Fe abundance of 
-0.29 ± 0.15 for 13 stars of their sample in NGC2004, in good 
agreement with the depletions of O, Mg and Si. Studies using 
F supergiants and Cepheids, however, report iron abundances 
representative of the mean metallicity (Hil l et al.|[1995[ [Luck[ 
et al.|1998| l. In this study, we adopt A[Z/H] ~ -0.3 ± 0.1 for 
the LMC. 

We note that for C and N the values in [Hunter et al.[P007| l 
are the best estimates of the baseline abundances in the Clouds. 
I.e. they are estimates of the C and N when the stars are born 
and before any rotational mixing has caused the photospheric 
abundances to change. It has been clear for some time that C 
and N are significantly underabundant in the Clouds in com- 
parison to the heavier elements. This may have an effect on the 
wind strengths of the hottest stars where CNO (mainly C) is a 
significant contributor to the line force, (see Table[T]i. 



4. Mass-loss determinations of early-type stars 

During recent years, three basic diagnostics to derive the mass- 
loss rate of early-type stars have been employed. These are: (a) 
infrared, millimetre and radio excess due to free-free processes; 
(b) ultraviolet resonance lines, and (c) optical lines, notably 
Ha, but for the hottest stars also Heii/i4686. For a didactic 
explanation of the way in which M can be derived from these 



diagnostics, see Lamers & Cassinelli ( 1999 1. In the context of 



The abundance of the a elements in the LMC relative to 
solar from |Hunter et al.[ ( [2007| l and [Trundle et al.| ( ,2007) is 



the dependence of mass loss on metal content, millimetre and 
radio excess emission in early-type stars is (at this moment) 
impracticable as the flux levels are weak and can only be mea- 
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sured for not too distant Galactic stars (up to a few kiloparsec), 
i.e. the method can not be appHed to Magellanic Cloud stars. 

UV resonance lines of relatively abundant elements such as 
C, N, O and Si are the most sensitive probes of mass loss, al- 
lowing detection of rates as low as ~ lO^^Mgyr"'. The lines 
typically saturate at about ~ 10"^ Moyr"^ therefore for strong 
winds they only provide lower limits to M. Because of the ease 
with which the lines saturate, unsaturated lines often relate to 
minor ionisation stages (and/or, obviously, to weak winds). The 
ionisation of these trace ions may depend rather critically on 
the treatment of line-blanketing, clumping, and shocks. Over 
the past decade, line-blanketing has been incorporated into 
model atmosphere codes for hot stars with winds. In a detailed 
comparison of three of these codes, |Puls et aT] ( |2005| l con- 
clude that the flux levels ai A > 400 A agree very well. Below 
400 A, discrepancies are found, implying that the population 
of species such as Niv, Civ, and Ovi should be taken with 
care, as should the mass-loss rates derived using these ions. 
One of the most important current challenges in model atmo- 
spheres treating outflows is to implement the direct and indirect 



effects of the line-driven instability (see e.g. Lucy & Solomon 
1970 Owocki 1994^ and references therein). This instability 



causes the formation of small scale density and velocity gradi- 
ents in the wind flow, creating "clumps" of gas and shocks lead- 
ing to X-ray emission and enhanced EUV-flux. We will discuss 



clumping in more detail in Sect. 8.2 Shocks may also impact 
on the ionisation of the species mentioned above. This presents 
a second reason for being cautious with the M values derived 
from UV resonance lines of which the ionisation is aff'ected by 
EUV and/or X-ray photons. A third reason to be cautious is that 
the abundances of C, N, and O may be affected by the surfacing 
of nuclear processed material (see Sect.|3]l. 

The third approach to derive the rate of mass loss is by fit- 
ting the Ho- profile, and, for stars of roughly spectral type 05 or 
earUer, He ii /14686. The advantages of fitting these lines com- 
pared to ultraviolet resonance lines are that /) the abundance 
determination is robust (for the most detailed approach, see 
Mokiem et al.||2065 2006 1; //) at least for hydrogen, the ion- 
isation balance is well known, i.e. it is not aff'ected by shock 
processes (which might still be a problem for helium). On the 
down side however, the Ha diagnostic is not as sensitive as are 
the UV lines. For OB stars, only for mass-loss rates in excess 
of ~ 10"^ Moyr"' can the contribution of wind emission to the 
Ha photospheric line be used to determine M. Even more, at 
such low wind densities the uncertainty about the velocity law 
(which cannot be constrained if the wind emission inside Ha 
is low) introduces an additional ambiguity, which increases the 
errors in the derived mass-loss rate considerably, up to factors 
of two to three (e.g., [Puis et al.|1996[ ). 

This poses a challenge: the Ha and UV resonance line di- 
agnostics have only a small mass-loss regime in common in 
which both techniques can be applied simultaneously. This 
is unfortunate in view of resolving the weak wind problem 
(see Sect. 8.1 1, although there may be alternative wind diag- 



nostics. First, the Bra line at 4.05 //m is intrinsically stronger 
than Ha and could push the sensitivity of the hydrogen line 



(as resulting from test calculations by JP and FN), providing 
sufficient overlap to derive M from both methods for a reason- 
able number of stars. Second there is a Pv resonance doublet 
at /i/illl8,1128 A. The abundance of phosphorus is about a 
factor of 10^ - 10^ less than that of C, N, O, and Si. As a re- 
sult of this the line does not saturate as easily. Consequently, 
it can be applied to stars with M in excess of 10"^ Moyr"' 
(see e.g. [Crowflier et"iL][2002l [Hillier et alT][2003] l. Moreover, 
for a range of O subtypes P v is expected to be the dominant 
ionisation stage (for recent results, see Puis et aLj2007) , mak- 
ing it less susceptible to shocks. Access to the far-ultraviolet 
region of the spectrum is provided by, for instance, the Far 
Ultraviolet Spectroscopic Explorer. Analysis of spectra ob- 
tained with this instrument strongly points to discrepant P v and 



Ha based mass-loss rates, suggestive of clumped winds ( |Massa 
let al.|2003t|Fullerton et al.|2006l l. 

5. Comparing M in different environments 

How should we compare the mass-loss rates of OB stars in dif- 
ferent galaxies? The most straightforward way of doing this 
would be to consider two stars with almost identical parameters 
- one in each galaxy - and directly compare their M values. 
There are two arguments against such an approach. The first 
one is of a practicable nature. The number of O and early-B 
stars studied in the Galaxy and the Magellanic Clouds is so far 
too limited to identify a significant number of stars with "iden- 
tical" luminosity, temperature, mass, and to a lesser extent rota- 
tional and terminal wind velocity. The impact of stellar rotation 
on mass loss appears important only for stars close to the Q-F 
limit (e.g. Maeder & Meynet 2000, and |Lamers|2004 and, e.g., 
Ivan Boekel et al.ll2003l and lSmith et al. 112003 1 for"aie extreme 
case of 77 Carinae), therefore, for relatively "normal" early-type 
stars one may consider it less critical. Terminal velocity diff'er- 
ences due to metallicity effects are found to be minor, both from 
observational (e.g. Evans et al. 2004b') and theoretical (e.g. Puis 



let al.|2000[|Krticka,2006) considerations. 

Even if for a few cases one could perform such a direct 
comparison, it would still not be an appealing approach, as the 
M(Z) relation may not necessarily be universal for all stellar 
parameters, and the result might not be applicable to all early- 
type stars. Therefore, the second reason against an object-to- 
object comparison is that it may ignore potential physical argu- 
ments that would allow the entire parameter space of hot mas- 
sive stars be used to establish the mass loss-metallicity relation, 
and that would add predictive power to the derived M{Z). 

Indeed the radiation-driven wind theory makes such pre- 
dictions, and a powerful way to proceed is through the use of 
the so-called modified wind momentum - luminosity relation, 
(WLR; e.g. |Kudritzki et al.|1995[|Kudritzki & Puls|2000 1 

logDmom = \og[Mvoo yfR) - xlogiUILo) + logDo . (1) 

In this relation the slope x and the constant Do are expected to 
vary as a function of spectral type and metal content (Kudritzki] 



et al. 1 1 999 1 1 Vink et al. 2000^ Puis et al. ,2000) . This equation ex- 



M-diagnostic by possibly a factor two to three ( Schaerer et al 



|1996) |Lenorzer et al.|2004( [Repolust et al.|2005 1 or even more 



presses that the mechanical momentum of the stellar wind flow 
is primarily a function of photon momentum. It is perhaps sur- 
prising that the stellar mass does not feature in this dependence. 
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The reason is that log Do contains a term oc (3/2 - x)logM, 
which (almost) vanishes since x happens to be ~ 3/2 for O- 
stars and early B-supergiants ( [Puis et al.|2000[ ). The uniqueness 
of the modified wind momentum relation has been confirmed 
by independent investigations (e.g |Vink et al.|2000}|Puls et al. 
2003 L In Tab. |3] coefficients for the WLR as predicted by Vink 
eraL|( [2000||2001) are given for the Galactic, LMC and SMC 



metallicities. 

Assuming the mass loss and terminal velocity are power 
laws of metallicity, i.e. 



M oc Z" 



and 



Zn 

it follows that 

(m + n) = AlogDmom/AlogZ 



(2) 



(3) 



(4) 



The index m may now be derived from a comparison of the 
modified wind momentum for different galaxies, adopting for 
instance the theoretical result n = 0.13 from ILeitherer et aL] 
1992| l to describe the Woo(Z) dependence. As the slope a of the 



WLR is not identical (though similar) for different metallicities 
one can not simply substitute the constant Do{Z) for Dmom in 
the above equation. What is required is a comparison of the 
actual Z),nom 3.1 a chosen luminosity, using the uncertainty in 



the WML relation at that specific L* (see Sect. 7.2 1. 



6. Observed mass-loss relations 

To determine the empirical mass loss versus metallicity de- 
pendence, we compiled mass-loss rates and terminal veloc- 
ity determinations from the literature, limiting ourselves to 
results obtained with state-of-the-art modeling techniques us- 
ing unified non-LTE line-blanketed stellar atmosphere models. 
Initially, we collected results obtained with four such codes: 
FASTwiND of |Puls et aT] ( |2005| l, cmfgen of |Hilher & Miller| 
('1998V WM-BASic of Pauldr ach et al.| ( |2001[ l, and isa-wind of 
[de Koter et al.| ( [1993, 199]}. We decided not to use studies 
performed with the latter two computer programmes (e.g. [de 



IKoter et al.''1994| [T998i [Bianchi & Garcia|[2002i [Garcia & 



|Bianchi^2004j ) as these account for only a relatively minor frac- 
tion of the total number of stars investigated (some 20 percent) 
of which some have been reanalysed with either fastwind or 
CMFGEN. This approach assures an extensive and relatively ho- 
mogeneous dataset, but not to the limit that only one code is 
used. In this way we can still investigate potential (systematic) 
differences between at least two methods. Objects for which 
the fastwind and cmfgen studies provided only upper limits on 
the mass-loss rate are not included in the determination of the 
empirical WLRs, however they are included in the figures for 
reference and comparison purposes. We also limit the sample 
to stars with T^ff > 24 kK. For lower effective temperatures, 
a decrease in the terminal wind velocity is observed (Lamers 



et al.|1995 Crowther et al.|2006] l, which may be accompanied 



by an increase in the mass-loss rate due to a change in the ion- 
ization balance of iron ( |Vink et al.||1999j , potentially leading 



to a different WLR for such relatively cool stars. This excludes 
some targets from the studies by Trundle et al. ( |2004 1; Trundle 
|& LennonI ( [2005[ ); [Crowther et al.| ( |2006| l. Finally, w e did not 
include the entries from the comprehensive study of |Markova 



et al. ( 2004) because, with the exception of the mass-loss rates. 



the stellar parameters were based on calibrations. We note how- 
ever that the stellar parameters derived by these authors com- 
pare well with other studies. 



6.1. Galaxy 

For the Galaxy we consider a sample based on the analyses 
performed by 



Repolust et al. (2004), Mokiem et al. (20051 



[Martins et"aL[ ( [2005| l and |Crowther et al., ( ,200 6). Note that the 



second study includes a reanalysis of the Cyg OB2 stars anal- 
ysed by'Herrero e"taL[([2002|, and HD 15629 and ( Oph studied 



by |Repolust et al. The relevant atmospheric parameters adopt- 



ing non-clumped mass-loss rates are listed in Tab. A.l Two 
entries are given for HD 15629 and HD 93250, as they were 
analysed separately using fastwind and cmfgen. In the follow- 
ing determination of the Galactic WLR, we adopt the results 
of the fastwind studies for reasons of consistency, but note that 
the differences of both analyses are not signficant. 

In Fig. [T] the distribution of the Galactic stars in the mod- 
ified wind-momentum vs. luminosity diagram is presented, 
where black/grey symbols refer to fastwind/cmfgen analyses. 
Different luminosity classes are distinguished using circles, tri- 
angles and squares for, respectively, class V, III-II and I objects. 
The open symbols show Dmom values resulting from scaling the 
M by a factor of 0.44, corresponding to -0.37 in log Dmom, for 
supergiants exhibiting Ha emission line profiles. A reduction 



of this amount was proposed by Repolust et al. ( 2004) to cor 



rect for the fact that these stars have a systematically higher 
wind momentum compared to O dwarfs and theoretical pre- 
dictions (see also 'P uis et"aL[[2003) parkova et alT[[2004| ). The 
physical interpretation for this systematic offset proposed by 
these authors is connected to wind clumping. Because Ha is a 
recombination line, its strength scales with the (wind)density 
p squared. In a uniformly clumped wind, the emission will in- 
crease by a factor / =< > / < p >^, where / is referred 
to as the clumping factor. As the Ha line-forming region of 
stars with Ha in emission is more extended compared to stars 
in which the profile is seen in absorption, these more extended 
regions must be more clumped than the innermost wind regions 
of those stars with absorption type profiles. In this interpreta- 
tion, therefore, the observed offset suggests a (spatial) gradient 
in the clumping factor or a difference in the clumping proper- 
ties of thin and thick winds. The (differential) clumping factor 
that corresponds to the applied scaling is / = 1 /0.44^ ^ 5. We 



will return to the issue of clumping in Sect. 8.2 



The left-hand side of Fig. [T] comp ares the stars that have 
been analysed by Mokiem et al. ( 2005 [ l in a homogeneous way 
by using an automated fitting method. We have constructed a 
WLR by fitting a power law, while accounting for both the sym- 
metric errors in L* and the asymmetric errors in Z)„iom, to the 
observed distribution. This empirical WLR is shown as a solid 
and dotted black fine for uncorrected and clumping corrected 
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Table 3. Coefficients describing empirical modified-wind momentum relations. Slope and vertical offset are given as x and Do, 
respectively. Parameters denoted with an apostrophe are derived from the wind momentum distributions including clumping 
corrected ZJniom values for stars with Ha in emission. For comparison, the coefficients of the theoretically predicted relation of 
Vinket al.| ( pOOOl l is also given. 



Galaxy sample 



log -Do 



MWG 



LMC 



SMC 



Mokiem et al. 
lotal 



2005 



Martins et al. 


20051 


Vink et al. 


20 


00 




Mokiem et al. 


2 


J07 



Mokiem et al. 
lotal 



2006 



1.86 ± 


0.20 


18 


71 ± 


1.16 


1.58 


±0.19 


20.16 ± 


1.11 


1.84 ± 


0.17 


18 


87 ± 


0.98 


1.56 


±0.16 


20.23 ± 


0.91 


3.15 ± 


0.95 


10 


29 ± 


5.08 










1.826 ± 


0.044 


18 


68 + 


0.26 










1.87 ± 


0.19 


18 


30 ± 


1.04 


1.49 


±0.18 


20.40 ± 


1.00 


1.96 ± 


0.16 


17 


88 ± 


0.91 


1.57 


±0.15 


20.02 ± 


0.84 


2.00 ± 


0.27 


17 


31 ± 


1.52 


1.50 


±0.23 


20.03 ± 


1.32 


1.84 ± 


0.19 


18 


20 ± 


1.09 


1.62 


±0.19 


19.26 ± 


1.10 




5.0 

log(L/Lsun/ 




5.0 5.5 
log(L/L3,,) 



Fig. 1. Modified wind-momentum - luminosity distribution of Galactic stars. Different luminosity classes are shown using circles, 
triangles and squares for, respectively, V, III-II and I class objects. Upper limits are shown as inverted triangles. Black/grey 
symbols refer to fastwind/cmfgen analyses. Clumping corrected Dmom values, using a clumping correction of -0.37 dex for 
objects with a Ha emission line profile, are shown with open symbols. Left: stars analysed by Mokiem et al. ( 2005[ l using an 



automated fitting method. The solid line corresponds to the empirical wind-momentum luminosity relation (WLR) for rates not 
corrected for clumping. A dotted line is the WLR for clumping corrected values. Right: total Galactic sample. Black solid and 
dotted lines, respectively, correspond to WLRs fitted to the complete sample using uncorrected and clumping corrected rates. 



The grey dashed line is the WLR obtained from fitting only the dwarfs analysed by Martins et al. ( 2005 1, adopting their rates for 
/ = 1 (grey circles). Upper limits were not considered in determining any of the best fits. 



wind momenta, respectively. The clumping corrected relation 
is found to be flatter because the clumping corrections only af- 
fected the two brightest objects. In Sect.[7]we will compare the 
empirical and theoretical WLRs in the observed range. 

On the right-hand side of Fig. [T] the observed Dmom distri- 
bution for the complete Galactic sample is shown. The empir- 
ical WLRs determined for this sample are again shown as a 
black solid and dotted line for the uncorrected and dumping- 
corrected rates. As can be seen in Tab.|3]the fit coefficients for 



both relations are in very good agreement with those derived 
from ffie homogeneously analysed sample. 

The Galactic sample comprises results from both cmfgen 
and FASTWiND studies. We note that fastwind is so far limited to 
studies of the visual spectral region - therefore He is the most 
important M diagnostics - while cmfgen (also) uses the ultravi- 
olet regime. In the latter approach the UV lines are given more 
weight in the mass-loss determination. In the case of weak 
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winds (< 10 ^Moyr ') the M determinations rely almost ex- 
clusively on fits to C IV ^^1548, 1551 (jMa rtins et al.|2005| l. 

Figure [T] shows good agreement between cmfgen and fast- 
wind studies for relatively high luminosities (log L^/L© > 5.5), 
consequently high wind densities. For lower luminosities, the 



UV analyses of the dwarf sample studied by [Martins et al. 



( 2005 | l show a systematic discrepancy with the average rela- 
tions. This is emphasised by the grey dashed line, which shows 
the average relation for these objects. The fit coefficients of 
this latter relation are also listed in Tab. [3] and clearly signal a 
discrepency between UV based mass-loss determinations and 
theoretical expectations. This is referred to as the "weak wind 
problem" (see e.g. Bouret et al.|2003| l. We will discuss this fur- 
ther in Sect.lOI 

6.2. LMC 



The studies from which we have drawn the LfVlC sample are 
those of [Crowther et al .1 (120021 , [Massey et al.| ( p004l p005) , 
[Evans et al.| ( |2004a| l and Mokiei neral.[ ( [2007| l. In Tab.|A.2[the 
relevant atmospheric parameters are listed. Note that some of 
the objects in this table have been analysed in more than one 
study. In those cases we adopted the results from studies using 



the automated fitting procedure developed by [Mokiem et al.| 



( [2005 j . If these were not available we preferred results from 
studies that used both the optical and UV spectral range over 
results that consider only the optical regime. For a number of 
stars analysed by Mokiem et al. ( 2007| l, no UV measurement 
of the wind velocity was available. Consequently, these authors 
estimated Uca by scaling the escape velocity at the stellar surface 
(fesc) with a constant factor of 2.6 (cf. [Lamers et al.|1995[ ). To 
account for the metallicity dependence and to facilitate a com- 
parison with theoretical predictions, we accordingly rescaled 
this Voa using Eq.[3]with Z = 0.5 Z© and n = 0.13 (cf. Leitherer 



et al.|1992) . In Tab. [A.2| the rescaled values are given between 
brackets. As also influences the density in the line forming 
region of wind sensitive lines (because of the requkement of 
mass continuity), also a rescaling of M was required. For this 
we used the wind-strength parameter 



Q 



M 



(5) 



that conserves the Ha equivalent width ( Schmutz et al.[[1989 
see also [Puls et al.[[T996| and [de Koter et al.[[1997| l. The com- 
bined effect of these re-scalings is a reduction of the modified- 
wind momentum by 0.08 dex. 

The distribution of the modified-wind momenta as a func- 
tion of stellar luminosity for the LMC stars is shown in Fig. |2] 
using the same symbols as in Fig. [T] On the left-hand side 
of the figure we (again) only consider objects that have been 
analysed by Mokiem et al. (2007) using an automated fitting 
method. The solid and dotted lines give the mean relations for 
uncorrected and clumping corrected rates. The correction ap- 
plied was the same as for the Galactic stars. We note however 
that the metallicity dependence of this correction is as yet un- 
known. 

For the total sample, shown on the right-hand side of the 
Fig. [2j the scatter is slightly larger In particular the bright su- 



pergiants Sk -67 22 at log L^/L© ^ 5.7 seems to stand out. 
Probably its discrepant position can be explained by the fact 
that [Massey et al.[ P005| could only determine a lower limit 
for its effective temperature, hence it could be intrinsically 
brighter. Despite the increased scatter, the obtained WLRs are 
in good agreement with the relations determined from the ho- 
mogeneously analysed sample. This can be seen from the fit 
coefficients given in Tab. |3] For a confrontation of the empiri- 
cal behaviour of the WLR with predictions, we refer the reader 
to Sect. [7] 



6.3. SMC 



In Tab. A. 3 atmospheric parameters are given for the SMC sam- 



pie that was compiled from the studies of [HilUer et al. 


20031), 


Bouret et al. ( 2003|l, Trundle et al. ( 2004|i, 


Massey et al. 


2004|, 


Evans et al.' ('2004a), 'Trundle & Lennon ( 


2005 1, Massey et al. 



in multiple studies we adopted results of one of these investi- 
gations following the same rules as applied to the LMC stars. 
Wind velocities given between brackets correspond to values 
that were calculated from the escape velocity at the stellar sur- 
face. These velocities and the associated mass-loss rates were 
scaled in a similar manner as was done for the LMC. Adopting 
Z - 0.2 Zq, the Dmom values were scaled down by 0. 18 dex. 

Figure |3] shows the distribution of the SMC stars in the 
^^mom VS. Li, diagram. The symbols used are the same as in 
Fig. [T] On the left-hand side of the figure only the objects 
that have been analysed using an automated fitting method by 
Mokiem et al. ( 2006) are shown. For log L^/L© < 5.3 the ma- 



jority of the mass-loss determinations are upper limits. The em- 
pirical WLRs are shown using a solid and dotted line for un- 
corrected and clumping corrected wind momenta. Note the po- 
sition of the dwarf NGC346-033 at log L*/L© ^ 5.0. We did 
not include this object in the fits, as its high wind momentum is 
the result of an anomalously high wind velocity resulting from 
the scaling with Uesc (also see |Mokiem et al.|2006| l. 

The right-hand side of Fig. [3] shows the wind momentum 
distribution for the total SMC sample. For this large sample 
the low luminosity part of the diagram still remains scarcely 
populated. Also note the upper limits determined using cmfgen 
for log L*/L© < 5.2. This could point to a steeper WLR relation 
for UV based mass-loss estimates compared to those obtained 
from Ha, as was found for the galactic case. However, this is 
not a firm statement as by far the bulk of the SMC targets in 
this luminosity range only show upper limits. 

The empirical WLRs for the total sample are shown on the 
right-hand side Fig.|3]as a black solid line for uncorrected wind 
momenta and as a dotted line for clumping corrected values. 
We find that these relations are strongly influenced by the ob- 
jects with log L*/L© < 5.3. This is shown by the grey solid and 
dashed lines, which correspond to the respective WLRs calcu- 
lated ignoring these objects. As the low luminosity part of the 
SMC diagram is rather uncertain, we opted to use this latter set 
of relations. Note that this choice will not influence our deter- 
mination of the metallicity dependence of M, as this wiU be 
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Fig. 2. Modified wind momentum - luminosity distribution for LMC stars. Symbols have the same meaning as in Fig.[l] Left: 
stars analysed by |Mokiem et al. { 2001) using an automated fitting method. Dotted and dashed lines are empirical WLRs fitted to 
uncorrected and clumping corrected ZJmom values, respectively. Right: total LMC sample. Fitted empirical WLRs are shown as a 
solid and dashed line for uncorrected and clumping corrected wind momentum rates. The fit coefficients for the different WLR 
relations are given in Tab.|3] Upper limits were not considered in determining any of the best fits. 




4.5 5.0 5.5 6.0 4.5 5.0 5.5 6.0 



log(L/L3,J log(L/L3,,) 

Fig. 3. Modified wind momentum - luminosity distribution for SMC stars. The symbols have the same meaning as in Fig. [l] To 
construct the uncorrected and clumping corrected empirical WLR (solid and dashed line) the dwarf NGC346 033 at log L^/L© ^ 
5.0 was not taken into account as the scaling of the escape velocity resulted in an anomalously high wind velocity. In Tab. [3]the 
fit coefficients for the individual relations are listed. Left: stars analysed by Mo kiem et al.| ( |2006| l an automated fitting method. 
Right: total SMC sample. It turns out that the WLR is strongly influenced by the objects with logL^/Lg < 5.3. This is shown by 
the grey lines, which correspond to the empirical WLRs calculated ignoring these objects. For establishing the empirical M(Z) 
as well as for the comparison with theory we will use the latter set of WLRs. Upper limits were not considered in determining 
any of the best fits. 



based on the stars brighter than L*/Lq= 5.3. In this range, the 
fit uncertainties are small (see Sect.[7|. 



7. Mass loss versus metallicity 

Now that we have established the empirical modified-wind mo- 
mentum luminosity relation for the three galaxies, we can de- 
termine the M(Z) relation. Before doing so, we first inter com- 
pare the results for these three environments, and confront them 
with predictions. 



7.1. Global comparison 

In Fig.|4]the empirical modified-wind momentum relations de- 
termined for the total observed samples (solid lines) are shown 



alongside the predicted relations of Vink et al. (2000 2001 



(dotted lines). The top, middle and bottom lines of each line 
style, respectively, correspond to the Galactic, LMC and SMC 
observed and predicted relations. To facilitate a meaningful 
comparison one sigma confidence intervals for the observed 
WLRs are shown as grey areas. First focusing on the left-hand 
side, showing the relations without a clumping correction, we 
see that the empirical relations are clearly separated from each 
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Fig. 4. Comparison of the observed wind momentum - luminosity relations (solid lines) with the predicted relations of 
et al. (2000 200 l) l (dotted lines). Top, middle and bottom lines of each line style, respectively, correspond to Galactic, LMC 
and SMC observed and predicted WLRs. Left and Right plots show, respectively, empirical WLRs including uncorrected and 



clumping corrected wind momentum rates. See Sect. 6.1 for the applied clumping correction. One sigma confidence intervals for 
the empirical relations are shown as grey areas. 



other beyond the fitting uncertainties. We interpret this as quan- 
titative evidence for a successive decrease of the mass-loss 
rates of massive stars in the Galactic, LMC, and SMC envi- 
ronment. 

For the entire investigated luminosity range, the relative 
separations of the empirical WLRs agree well with the sepa- 
rations predicted by Vink et aL] ( [2001) , although the empirical 
WLRs do show a systematic offset compared to the theoreti- 
cal results. Focusing for a moment on the relative separations 
only: measured at logL^/L© - 5.75, which coincides with the 
region where the fit uncertainties are the smallest, the offsets 
between the empirical and theoretical relations relative to the 
Galaxy are: 0.28 and 0.25 for the LMC, and 0.62 and 0.57 for 
the SMC. 



The systematic offset between observations and theory is 
of the order of 0.2 dex. This suggests that ;/ wind clumping 
does affect the determination of empirical mass-loss rates, but 
would only have a marginal effect on mass-loss predictions , the 
empirical mass-loss rates would be overestimated by at most a 
factor of ~ two. 



Indeed, when turning to the clumping corrected empirical 
relations, the relative offset almost disappears. This is shown 
on the right-hand side of Fig. |4j where we show the wind mo- 
mentum relations obtained from the total observed samples ac- 
counting for this correction. Note that, strictly speaking, this 
agreement pertains only to stars not showing Ha in emission, as 
the clumping corrections is designed such as to scale the ones 
that do show Ha emission to those that do not (see Sect. |6.1| l. 
The slopes of the empirical relations, however, are slightly flat- 
ter, as the clumping correction preferentially affects the high- 
luminosity objects. 



7.2. The empirical MiZ) relation 

To calculate the mass loss metallicity dependence we use the 
relative separation between the empirical modified-wind mo- 
mentum relations at log L^^/Lq = 5.75, for which the fit un- 
certainty in all three relations is at its minimum. Note that 
the differential slopes also allow for a luminosity dependent 
Z-dependence, however, based on the uncertainties we doubt 
whether this would be meaningful. Moreover, for the SMC it 
was shown that this slope is very sensitive to the rather un- 
certain low luminosity domain. In Tab.[3]the power law indices 
for the observed M(Z) dependence for the LMC and SMC rela- 
tive to the Galaxy are listed. These were calculated using Eq.|4] 
assuming n = 0.13. The indices determined from the homo- 
geneously analysed sample are in good agreement with those 
determined for the total sample. For the LMC relatively larger 
differences are found compared to the SMC, which is a result 
of the smaller metallicity difference of this system relative to 
our galaxy. This also explains the larger error bars. 

To determine the global metallicity dependence, we made 
a linear fit regarding the logarithm of the differential (i.e. with 
respect to the Galaxy) modified wind momentum and the log of 
the differential metallicity. Assuming smooth winds, we find: 



MocZ 



0.83±0.16 



(6) 



We emphasize that the error bar in m is a fitting error of 
the mass-loss rate and the metallicity only. Several uncertain- 
ties add to the quoted error These include /) those related to 
the adopted luminosity at which we gauge the relation; //) er- 
rors related to the neglect of the correlation between luminosity 
and modified wind momentum ( |Markova et al.|2004l see); 
errors related to uncertainties in the Vao{M) dependence, i.e. in 
the adopted value n - 0.13, and iv) those related to possible 
systematic uncertainties in the metallicities of the Magellanic 
Clouds. The motivation for gauging at log L^/Lq - 5.75 is 
given above. Least squares fit analysis including the correla- 
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tion between L and D yield similar results for the stars anal- 
ysed using the automated fitting method, though for the total 
sample the Galactic log D{\og L) relation is found to be some- 
what steeper. Deviations in the Galactic slope are to be ex- 
pected ( [Markova et al.||2004| l. However, at log L,^/Lq = 5.75 
the wind momenta are not found to be very different and result 
in uncertainties Am < 0.10 dex. Uncertainties in the velocity 
vs. metallicity dependence are not very well known, however, 
they are likely to be small as well. By far the largest source of 
error in m is that due to uncertainties in the metallicities of the 
Magellanic Clouds, notably that of the SMC. To illustrate this, 

if the SMC metal content is not A[Z/H] 0.7 + 0.1 but rather 

-0.5 + 0.1, the empirical value of m will increase to 1.13 + 0.25. 

Keeping this in mind, the derived exponent m of the empiri- 
cal mass-loss luminosity relation as given in Eq.|6]appears to be 
consistent with the predicted value of mpied = 0.69 + 0.10 ( Vink 
et al.|2001)) and /Wpred = 0.67 from the alternative investigation 



by Krticka 
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If we apply the earlier-mentioned clumping correction scal- 
ing of -0.37 dex in logDmom for stars showing Ho' emission, 
the metallicity dependence is found to be: 



MocZ' 



■0.72±0.15 



(7) 



This result too is consistent with both theoretical predic- 
tions quoted above. 

Three remarks need to be made. First, we can only state 
that this M(Z) relation holds for stars more luminous than 
~10^'^ Lq. The reasons are the lack of firm M determinations 
for stars less luminous than this number and the inconsistency 
between the Ha and UV wind lines mass-loss diagnostics (see 
below). Note in particular that the slope of the WLR might de- 
pend on the wind-density and/or metallicity itself (as it is pre- 
dicted from a close inspection of the hne-strength distribution 
functions, cf. |Puls et al.|2000 1, though not to such an extent as 
indicated by the weak wind problem (cf Fig. [T]). In so far, it 
is possible that the metallicity dependence for low-luminosity 
objects deviates from our previous results. 

Second, the absolute agreement in wind momentum only 
holds if the winds of the O dwarfs are not significantly 
clumped. If they are, the predictions (which do not yet ac- 
count for clumping) may overestimate the mass-loss rate and 
therefore Dmom- Third, ifO dwarf winds are clumped, but the 
clumping does not depend on metal content or mass-loss (i.e. 
/ f{Z, M{Z))) then the derived m will not be affected, and, 
therefore, the agreement in the observed and predicted mass 
loss - metallicity scaling will remain preserved. To rephrase 
this last point: irrespective of what the actual M values are, if 
clumping were independent of metallicity and wind density, the 
derived power-law index is anticipated to be correct - even if 
the winds are significantly clumped. 

Interestingly, these analyses suggest that the systematic er- 
rors in the slope of the mass-loss-metallicity relation, m, as in- 
troduced by potential wind clumping may be less relevant than 
remaining uncertainties in the metallicity determinations of the 
Magellanic Clouds. 



8. Discussion 

8.1. The weak wind problem 

As mentioned before, for the majority of stars at lower lumi- 
nosities the ultraviolet resonance lines based mass-loss rates 
disagree strongly with theoretical predictions. Still, there are 
few (though different) Galactic objects at similar, low luminos- 
ity which display a "normal" Ho- mass-loss rate, which might 
point either to different physics or questions the validity of at 
least one of the applied diagnostical tools. 

Going from high to low luminosity, figures [T] and [3] show 
that this discrepancy appears at » 10^^^ Lq. The dwarf star 



sample recently analysed by Martins et al. ( 2005 i using cmf- 
GEN suggests progressively weaker winds than are predicted by 
current theory. At L ~ 10"^^^ Lq this discrepancy has reached a 
factor 100. Note that given the still limited number of stars in- 
vestigated at relatively low luminosity alternative descriptions 
(compared to a steeper slope) of the WLR might also be ap- 
propriate. For instance, the WLR may retain the same slope 
but jump to a two order of magnitudes lower value of Dmom 
at L < 10^^^ Lq compared to the predictions. In the SMC, the 
discrepancy could be even worse, as UV analyses only yield 
upper limits. Note that also the Ha analyses of SMC objects 
at L < 10^^^ Lq yield almost exclusively results as upper lim- 
its. So, in principle, these need not be contradictory to the UV 
based M. 

What could be the reason for this "weak wind problem"? 
There are three directions in which one might look for causes: 
first, errors in the Ha and/or UV analysis; second, missing 
physics or invalid assumptions in the wind predictions, or, 
third, the nature of the stars showing the weak wind problem 
is different from that of normal OB stars (that do appear to fol- 
low theory). Here, we focus on the first possibility. For more 



thorough discussions on this topic we refer to Martins et al. 
( |2004||2005| l and |de Koter| ( [2006| . 

Two important advantages of the Ha method are that the 
ionisation and abundance of hydrogen are well known. One 
cannot make this statement with similar confidence for the el- 
ements responsible for the UV resonance lines. For the weak 
wind stars, the prime M diagnostic - C iv as mentioned - is only 
a trace ion, i.e. it does not represent the dominant ionisation 
stage. This makes it is susceptible to X-ray emission, thought 
to originate from shocks that may develop in the wind outflow. 



Martins et al. ( 2005 i present test calculations estimating that 
X-rays may decrease the C iv ionisation by an order of magni- 
tude, requiring an increase of the mass loss by the same amount 
to preserve the line fit. With typical uncertainties in the carbon 
abundance of a factor of two, this shows that the UV method 
is much more prone to errors than the Ha method. This does 
not imply that Ha does not suffer from uncertainties. Besides 
the ambiguity in M introduced by the unknown velocity law 
in weaker winds which has been discussed already, continuum 
rectification issues and potential nebular emission limit the ap- 
plicability of the Ha diagnostic to stars with M > lO^^Moyr'^ 
To reach this limit requires high signal-to-noise data of stars 
with no appreciable nebular emission, such as is the case for 
three of the five Galactic stars at L < 10^ Lg - i.e. in the (UV 
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identified) weak wind regime. These galactic stars are ^Oph, 
CygOB2#2, and HD 217086. For two other stai's at these low 
luminosities, rSco and lOLac, we also established the mass- 
loss rate, however, the error bars on these determinations are 
large. For the three stars mentioned, the Ha derived rates do 
appear to be more or less consistent with predictions. Note that 
some of the C iv analyses account for (canonical) X-ray emis- 
sion suggesting that the problem may not be fully explicable 
by uncertainties in the UV method. 

So far the origin of the weak wind problem has not yet 
been identified. We stress that because of the poor overlap be- 
tween the Wa and UV diagnostics, it is at present difficult to 
exclude the possibility that the (major) cause is simply con- 
nected with uncertainties in the wind structure predicted by the 
model atmospheres that are needed to derive the empirical M- 
notably the ionisation of species responsible for the UV reso- 
nance lines. Further spectroscopy in the IR, particularly using 
Brc as an ideal mass-loss indicator also for weak winds, might 
unravel at least part of the problem. 

8.2. Small-scale clumping 

So far, we have not discussed the physical reality of the "clump- 
ing" correction applied in Sect. [6] to bring the WLR relation 
of the supergiants with Ha emission line profiles in accord 
with those of dwarfs. Is wind clumping observed in O stars? 
And, is clumping more important in the winds of supergiants? 



Clumping was studied in detail by HilUer ( 1991 1 to explain the 
scattering wings of He ii lines of Wolf-Rayet stars. Although 
these stars are known to have very dense winds, their He ii scat- 
tering wings are weak features. Consequently, for the weaker 
O star winds these diagnostics cannot be applied. Nonetheless, 
direct observational evidence that the winds of O stars are also 



clumped was provided by [Eversberg et al. ( 1998 1 from time 



series spectroscopy of He ii /14686 in f Pup, and by Crowther 

Bouret et al.|p003|);|Massa 



et al. 


(20021 


et al. 


(|2003]l 



HiUier etal.|(|2003]l; 



Martins et al. 



( |2004l 



2005] l; |BOTget et al.| ( |2005] l. 



and FuUerton et al. (2006 ) from analysis of UV (resonance) 
lines, partly in combination with optical lines. 

The fitting of the ultraviolet lines, notably 
Pv/l/11118, 1128, Ov/11371, and NiV/11718, is found to 
be improved if a distance dependent clumping is introduced 
- with clumps starting to form above the sonic point and 
reaching maximum clumping in the outer winc|^- and a void 
interclump medium is assumed. In the above listed studies a 
total of 25 stars have been analysed using the cmfgen code. 
Clumping corrections are applied in roughly half the sample, 
including both dwarfs and supergiants - for objects brighter 
than IO^ '^Lq. The maximum clumping factors that are derived 
from these studies are in the range 10 - 100. Though clumping 
corrections have turned out to be necessary for some stars 
but not for others, this does not exclude the possibility that 
all stars are clumped to a certain degree: Due to the fact that 
only few lines react differentially on clumping and thus allow 



' In this respect "outer" implies the formation region of the UV 
resonance lines, which is inside of the region that contributes the bulk 
of the radio flux. 



for distinct statements, it is rather difficult to exclude the 
presence of clumping in a certain object at the present state of 
knowledge unless the wavelength coverage is fairly complete 
(extending from the UV to the IR), and unless a large number 
of systematic investigations have been performed. 

In addition to these purely spectroscopic investigations. 
Puis et al. ( 2006| l analyzed a sample of Galactic O-type (su- 
per)giants more recently, combining Ha, infrared, mm and ra- 
dio fluxes to derive constraints on the radial stratification of the 
clumping factor. Since all diagnostics employed have a de- 
pendence, only relative clumping factors could be derived, nor- 
malized to the values in the outermost, radio-emitting region, 
whereas absolute values for the clumping factors and thus the 
actual mass-loss rates remained unconstrained. 

In contrast to present hydrodynamical simulations of the 
line-driving instability but in (partial) agreement with previous 
studies it was found that for denser winds clumping starts fairly 
close to the sonic point, reaches a certain maximum (mostly 
within r < 2Ri^) and decreases outwards, where the maximum, 
normalized clumping factors are of the order of 3 to 6. For 
weaker winds, on the other hand, the clumping factors in the 
inner and outermost regions turned out to be similar 

This appears to be consistent with findings from the linear 
polarimetry of Luminous Blue Variables ( [Davies et al.||2005] ), 
where it was found that (i) clumping must start close to the 
photosphere to reproduce the observed levels of polarimetry, 
and (ii) denser winds yield a larger amount of linear polarisa- 
tion (see also |Robert et al.|1989] for WR stars). 

The results by [Puis et aT] ( |2006[ ) could mean one of two 
things: either the outer wind is significantly clumped whilst the 
inner wind is even more strongly clumped (which would im- 
ply that the UV results calling for large clumping factors are 
correct), or the outer wind is not significantly clumped and the 
current mass-loss predictions are approximately coiTect. The 
latter possibility, of course, would require an explanation of the 
different UV results. 



8.3. Concluding remarks 

The first result we mention is that the empirical M(Z) relation 
appears to be in agreement with theoretical predictions for lu- 
minosities larger than ~10^^ L©. Second, for lower luminosi- 
ties, UV based M determinations for dwarfs seem to indicate a 
breakdown of the theory, leading to discrepancies of up to a fac- 
tor 10^ for the lowest luminosity stars with a measurable mass 
loss. We found that for the Hor based mass-loss rates of three 
galactic stars, ^Oph, CygOB2#2, and HD217086, do follow 
the predictions down to log L^/L© a; 4.9. Their M values are 
~ 2 X 10"^ Moyr"' or less, which is so low that they may suf- 
fer from uncertainties that are not reflected in the derived error 
bars, for instance those associated with the uncertain velocity 
law in low density winds and/or the unnoticed presence of mi- 
nor amounts of nebular emission. Still, given these findings, it 
appears premature to exclude the possibility that the (major) 
cause of the breakdown between empirical and predicted mass 
loss at low luminosity is connected to modeling uncertainties 
associated with the empirical M values. 
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To further probe the cause of the possible theoretical break- 
down at low luminosity, it is therefore necessary to study those 
stars for which both the Ha and UV line indicators of mass 
loss can be applied - although unfortunately there is only a very 
limited mass-loss regime for which this is possible. To this end, 
the use of Bra as a mass-loss diagnostic ( jLenorzer et al.|2004{ 
Repolust et al.|[2005j ) will prove to be fruitful as it will widen 



the overlap between hydrogen and UV line analyses. Regarding 
clumping, the Pv resonance line seems a powerful diagnostic 
( |Massa et al.|2003[[Fullerton et al.|2006[ ). Also, it may be fruit- 
ful to focus on the question whether all O-type stars suffer from 
clumping. Concerning the effect of clumping on the mass loss 
vs. metallicity relation, the following important remark can be 
made: regardless what the actual M values are, if clumping is 
universal in O-type stars, and if it does not depend on wind 
properties such as density or metallicity, then the derived scal- 
ing (i.e. the power law index m) remains correct. 

Finally, let us briefly contemplate some implications of O- 
type star mass-loss rates that are lower by factors of 3 to 100 
than so far assumed. These implications branch out into at least 
two directions: wind hydrodynamics and stellar evolution. For 
our understanding of the intrinsic instabilities associated with 
the wind driving mechanism (Owocki 1994\ the requirement 
of clumping factors / ~ 10 - 100 in the Ha line forming re- 
gions of dwarf O-type stars, which do not extend far beyond the 
sonic point, seems to imply that these density inhomogeneities 
have developed already in the sub-sonic part of the flow. As an 
example, the 05 V star Nl 1-051 in the LMC has an Ha line 
forming region that extends from the photosphere to a velocity 
of 13 kms"', which is about 0.6 times the sonic velocity. For 
other O-type dwarfs, similar results are found. This is not an- 
ticipated by theory, which predicts that the onset of line-driven 
instabilities is only at about the sonic velocity (where, admit- 
tedly, the growth rate of the instability is large). Implications 
for the wind driving mechanism may also be severe, although 
it should be noted that clumping has not yet been incorporated 
in mass-loss predictions. 

Clumping may also impact on our understanding of mas- 
sive star evolution. Weaker stellar winds will cause less loss of 
angular momentum. Consequently the stars will not spin down 
as rapidly as currently thought. It may even be expected that 
most Galactic massive stars retain their initial rotational ve- 



locity properties during the entire main sequence (Meynet & 



Maeder|2000}|Maeder & Meynet|200T) . For flie brightest stars 
- corresponding to initial masses of 60 Mq or more - the inte- 
grated main-sequence mass loss could drop dramatically. With 
the current rates they are expected to lose 20 to 40 percent of 
their mass in the H-burning phase. This could drop by a factor 
of three, requiring the objects to lose ten to tens of solar masses 
by other means in order to evolve towards the hot Wolf-Rayet 



phase (Smith & Owocki 2006). This could either imply an in- 
crease of the duration of the LBV phase by a factor of two, 
assuming that the (unknown) mechanism thought to be respon- 
sible for the eruptive mass loss in this phase is unaffected by 
clumping issues (see e.g. Humphreys & Davidson|1994) l, or al- 
ternatively that some LBVs explode before reaching the Wolf- 
Rayet phase ( |Kotak & Vink|2006| |Gal-Yam et al.|2007| |Smith| 
[2007] l. 
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Table A.l. Wind parameters of Galactic O- and early B-type stars. 
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Table A.2. Wind parameters of O- and early B-type stars in the LMC. Wind velocities given between brackets are calculated 
from the escape velocity at the stellar surface. Identifications: "BI" from Brunei et al. ( 1975 1, "LH" from Lucke ( 1972), ex cept 



LH 51-496, which is identified by Gamany et al. ( 1994), "Nil" from |Evans et al.| ( |2006| l,^136" from [Hunter et al.| ( |1997| l and 



Massey & Hunter| ( |1998| ) and "Sk" from ,5 anduleak, ( ,1970^ . 
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Table A.3. Wind parameters of O 

from the escape velocity at the stellar surface. Identifications: "AzV" from 



and early B-type stars in the SMC. Wind velocities given between brackets are calculated 

Azzopardi & Vigneau| ( [1975] [1982) 1, "MPG" from 



[Massey et al.| ( [T989| l and "NGC330" and "NGC346" from |Evans et"ai:] ( [2006 
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